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1 Relativistic quantum field theories for complex scalar and Dirac fields are investigated 

in non-equilibrium thermo field dynamics. The thermal vacuum is defined by the Bo- 
goliubov transformed creation and annihilation operators. Two independent Bogoliubov 
parameters are introduced for a charged field. Its difference naturally induces the chemi- 
cal potential. Time-dependent thermal Bogoliubov transformation generates the thermal 
counter terms. We fix the terms by the self-consistency renormalization condition. Eval- 
uating the thermal self-energy under the self-consistency renormalization condition, we 
derive the quantum Boltzmann equations for the relativistic fields. 
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1. Introduction 

Much attention has been paid for thermal evolution of high energy phenomena, for 
example, hot and dense quark matter. It is expected that the time dependence of 
the system can be observed in a relativistic heavy ion collisions and cooling process 
of dense stars. In such a system a thermal quantum field theory, i.e. a quantum field 
theory to describe thermal dynamics, is necessary to investigate the non-equilibrium 
phenomena for relativistic fields. There are several theories, the the rmo field dy- 
namics (TFD), the Kadanoff-Baym formalism, the Langevin equatioiffEmill and 
so on. In this paper we focus on thermo field dynamics to describe a non-equilibrium 
thermal system for relativistic charged fields. 

TFD has been proposed by Y. Takahashi and H. 

UmezawEPSEEl. in TFD the 

canonical formalism for a quantum field theory is established at finite temperature. 
The so-called thermal vacuum state is defined to calculate the thermal average 
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2 Y. Mizutani and T. Inagaki 

of dynamical operators. It introduces the thermal Bogoliubov transformation. The 
fundamental formalism of TFD for relativistic complex scalar and Dirac fields is 
developed in a finite and homogeneous system, see for example Ref. ]9\ A propaga- 
tor for the Dirac field has been discussed at finite chemical potential in Ref. [101 In 
Ref. [11] the complex scalar field has been applied to TFD in a spatially inhomoge- 
neous system, and the thermal conductivity and the diffusion coefficient has been 
calculated for interacting pions in hot and dense nuclear matter. 

TFD has been extended to describe an out of equilibrium systenl . In the 
non-equilibrium thermo field dynamics (NETFD) the thermal Bogoliubov trans- 
formation is generalized to depend on time variable. The time -depende nt thermal 
Bogoliubov transformation induces the thermal counter term | 13 | 14 | 15 | rpj^g term 
is fixed by imp osing the self-consistency renormalization condition at the equal 
time limil! 1 -^ ^. Then the quantum Boltzmann equation is derived from the renor- 
malization condition. This procedure is a pplied to in vestigate the Bose-Einstein 
condensation in trapped cold atom systems ^ 6 | 1 7 | 1 8 | l H nexFD is also extended to 
an inhomogeneous system with diffusion'^^ 

In Ref. [22] NETFD has been applied to a neutral scalar field based on the canon- 
ical quantization. In NETFD the neutral scalar field depends on a thermal Bogoli- 
ubov parameter, which is equal to uncharged particle number density under the 
self-consistency renormalization condition. The relativistic charged fields depend 
on two types of the Bogoliubov parameters, n^(t) and n^{t), which correspond to 
the particle and anti-particle number density under the self-consistency renormal- 
ization conditions, respectively. The thermal counter term for the neutral scalar field 
fixes one of the Bogoliubov parameters. An alternative condition or assumption is 
necessary to fix all the Bogoliubov parameters. In the case of the Dirac fermion it 
has been pointed out that the positive and the negative frequency parts can not be 
described in an usual Lagrangian form and conjectured that the charge conjugation 
is broken in the framework of NETFD I^Sl. 

In this paper we approach relativistic complex scalar and Dirac fields at finite 
chemical potential based on the canonical quantization. In §2 we consider a scalar 
field in NETFD. Following the canonical formalism, we quantize the complex scalar 
field with finite chemical potential. The time-dependent Bogoliubov transformation 
modifies the unperturbed Hamiltonian and introduces the thermal counter term. 
In §3 we discuss the self-consistency renormalization condition for the complex 
scalar field. The time evolution equations are derived under the condition. In §4 a 
two loop quantum correction is calculated in a {(f)^ (f))"^ interaction model. Applying 
the self-consistency renormalization condition at the lowest order, we find that the 
time evolution equations for the Bogoliubov parameters coincide with the quantum 
Boltzmann equations for the charged Bose particles. In §5 we study a Dirac field in 
NETFD. The Dirac field can be quantized at finite chemical potential following a 
similar procedure with the complex scalar field. Because of the time-dependence of 
the thermal Bogoliubov transformation, the unperturbed Hamiltonian is modified 
and the thermal counter term appears. In §6 we extend the self-consistency renor- 
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malization condition to the Dirac field. The time evolution equations for the Dirac 
field are obtained from the condition. In §7 we apply the result to a Yukawa inter- 
action model. One loop quantum corrections are calculated for the neutral scalar 
and the Dirac propagators. Imposing the sclf-consistcncy renormalization conditions 
on the scalar and Dirac fields, we reproduce the well-known quantum Boltzmann 
equations at the week coupling limit. Concluding remarks are given in §8. 

2. Canonical Quantization for a Complex Scalar Field 

We start from the description of a complex scalar field at finite chemical potential 
in NETFD. In this paper we confine ourselves to homogeneous processes. Thus the 
complex scalar field is expanded in terms of two independent operators, Up and bp, 
where p denotes the spatial momentum of a single Fourier mode. The creation and 
annihilation operators obey the commutation relations, 



(27r)35(3)(p-fc), (1) 
[bp,bi] ^{27rr S(^Hp~k), (2) 
[ap, bl] = [al bl] = [ap, bk] = Wp, bk] = 0. (3) 

In TFD the tilde operators, Up and bp, are introduced according to the tilde 
conjugation rules, 

{A,A2y^A,A2, (4) 
{ciAi+c2A2y=clAi+4A2, (5) 
(ir= A, (6) 
{A^r=A\ (7) 

where Ai and A2 show any operators, ci and C2 are c-numbers. Following the tilde 
conjugation rules, we find that the tilde operators, ap and bp, obey the commutation 
relations, 



[ap,al] 



(2^)3<5(3)(p-fe), (8) 

[bp,bi] = {27:)'5^'\p~k), (9) 

[ap, bl] = [al bl] = [ap, bk] = [a^ , bk] = 0. (10) 
The ordinary operators commute with the tilde operators, 

[ap, ttk] = [ttp, al] = [up, bl] = [up, bk] = 0, (11) 

[bp, bk] = [bp, bl] = [bp, ak] = [bp, al] = 0. (12) 

Time evolution of the field is described by the hat-Hamiltonian which is constructed 

by 

H = H-H, (13) 



where H is an ordinary Hamiltonian and H is the tilde conjugate Hamiltonian. 
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4 Y. Mizutani and T. Inagaki 

In TFD the thermal expectation value of a dynamical operator is given by the 
expectation value under the state so-called thermal vacuum state, {6\ and \9). For 
the complex scalar field the thermal vacuum state is defined by the two operators, 

and 77p, and their tilde conjugates, and rip, which are the thermal Bogoliubov 
transformations of Op, bp, hp and hp with two independent thermal Bogoliubov 
parameters, n±{t]p), 

^^{t) = B{n+{t-p)rf^al{t), W = -al{t)B-\n+{t-p)f-, (14) 
<(i) = B{n^{t-p)rHl{t), fj^it) = mt)B-\n,{t;p)r-. (15) 

Here and below the upper indices indicate the thermal doublet notation. In this 
notation creation and annihilation operators are described as 



hi 



~^p={al-ap)- (16) 



The Bogoliubov transformation has a 2 x 2 matrix form. 

In the non-equilibrium and homogeneous system the thermal Bogoliubov pa- 
rameters, n±(t;p), depend on time, t, and the magnitude of the momentum, p. We 
set the thermal Bogoliubov matrices, B and B~^, as 

Bin^)=('^:^~:^), (17) 



-1 1 



(18) 



This expression makes available to calculate quantum corrections in the Dyson- Wick 
formalisnP31. 

According to the transformed operators, £,p,rip,£,p and fjp, the thermal vacuum 
state is defined as 

^p\d) = ip\d) = Vp\0) = Vp\0) = 0, (19) 

(^1^^ = ml = {0\vi = ml - 0. (20) 

It should be noticed that the physical observables are constructed by the original 
operators, Qp and bp. 

In an equilibrium system the thermal vacuum state has a condensate nature 
structure with the thermal pairs, aa and bb, ^ 

1^) = UBAXapBMXbM), (21) 
(^1 = mUB!aiXa)Usl{Xb), (22) 

where |0)) and ((0| represent the non-thermal vacuum states, 

ap|0)) = ~ap\0)) = bp\0)) = bp\0)) = 0, (23) 
((0|at = ((0|5t ^ ((0|5t = mil = 0. (24) 
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The transformation operators are described as 

UbAXu) = exp[ixQiG'ai]exp[-iXa2G'a2]exp[jXa3G'ai]exp[zxQ2GQ2]exp[-iXaiGai], 

(25) 

UsfiiXb) = exp[ixfciG6i]exp[-ix62Gb2]exp[zxfc3G6i]exp[zxb2Gh2]exp[-ixfciG6i], 

(26) 

with the generators, 

Gai{k) = i{a\ak + a\ak), Ga2{k) = i{akak + a\al), (27) 
Gbi{k) = i(blbk + blbk), GUk) = i{bkbk + blbl). (28) 

We use an abbreviation, 

Xa{b)iGa{b)j = j d^k Xa(b)iik)Ga(b)j{k). (29) 

The thermal vacuum state is classified according to the parameters, 
Xai,Xa2,Xa3,Xbi,Xb2 and Xb3- From Eqs. and we can find that the ex- 
pectation value of operators under the thermal vacuum coincides with the ther- 
mal trace of the operators. By using the transformation operators, UB,a{Xa) and 
UB,b{Xb), the thermal Bogoliubov transformations are written as 

C = UBAXa)a;Us^,{xa), S.; = UBAXa)a'iUs\{xa), (30) 
r?^ = UB.biXbKU^lixb). = UB.b{XbWpUs\{xb). (31) 

The thermal Bogoliubov parameters, n+(fc) and n-{k), can be represented as a 
function oi Xai{k) and Xbi{k), respectively. However, the operators, UB,a and UB,b, 
contain divergent coefficients in field theories with an infinite degree of freedom. 
These operators are mathematically ill-defined. Here we employ well-defined ex- 
pressions p4)) and (fTS)) and quantize the fields. 

The Bogoliubov transformation keeps the commutation relations. Thus the 
transformed operators, and rjp^ and their tilde conjugates, and f)p, satisfy 
the commutation relations, 

K.,eI] = Kp,el] = (2^)^<5(-^nP-fe), (32) 
[VpA]^[%^vl] = {'2^f5^^\p-k), (33) 

with all other commutators vanish. The scalar field is quantized under these com- 
mutation relations. 

We develop the quantum field theory for the complex scalar field with finite 
chemical potential according to NETFD which has been proposed by H. Umezawa 
and Y. YamanakaP^. Below we work in the interaction picture. To guarantee the 
time independence of the thermal vacuum we assume that the complex scalar field 
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6 Y. Mizutani and T. Inagaki 



is decomposed as a Hermitian form in terms of the transformed operators, 

d^p 1 



(2^)3 
d^p 1 
(2^)3 72^ 



^ipx 



with 



and 



(34) 
(35) 



(36) 
(37) 

(38) 
(39) 



where ojp is the relativistic energy eigenvahie,\/p^ + m^, and fi denotes the chemical 
potentiaL The energy eigenvalue, ujp, generally depends on time through the screen- 
ing effects and so on. In this paper we focus on a system with a time independent 
energy eigenvalue, for simplicity. The chemical potential is introduced according to 
the conserved charge in an equilibrium system. See appendix A. 

The tilde conjugate field is defined by the tilde conjugation rules, (jl])-®. 

Thus the time dependent operators can be combined as the thermal doublet nota- 
tion, 



The complex scalar field is written in the thermal doublet notation. 



(42) 
(43) 



d^p 



1 



(27r)3 
d3p 



jpx 



f)p{tx) 



-ipX 



(27r)3 
d^p 



{^pit.)e'P-'' + {rs%{t.fre-^P-^}, 



(44) 



(27r)3 
d^p 



-ipX 



(27r)3 



{Cp{tx)e-'P-'^ + {Vp{txrr,re^P-^], 



(45) 
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Non- Equilibrium Thermo Field Dynamics for Relativistic Complex Scalar and Dirac Fields 7 

where t^^ is the third Pauli matrix acting on the thermal indices. The canonical 
conjugate field is given by 



The canonical commutation relations are satisfied at the equal time, 

[cj,^{t,x),7:l{t,y)] = [^^it,x),7r^^it,y)]=zS'^^\x-y)d'^^. (48) 

Next we decompose the complex scalar field in terms of the original operators, 
Up and bp. Differentiating the operators, ^p{tx) and £.p{tx), with respect to the time 
variable, we obtain 

dtj^it,) = -lu^+Apii-.)^ (49) 
da'^{tx)=iLo+J'^{tx). (50) 

The thermal Bogoliubov transformation (|14p gives the time evolution equations for 
the operators, ap{tx) and a'^{tx), 

dt.a^itx) = -^W+.p - in+(t.;p)To}"^a^(t,), (51) 
dt^a^it,) = a^(t,)z{w+,p - zn+(t,;p)ro}'^", (52) 



where the matrix, Tq^, is given by 



nap / 1 -1 



To'-[l_l)- (53) 

It should be noted that the matrix, T^^, satisfies, T^^T^'' = 0. Eqs. 1^ and 1^ 
indicate that the energy eigenvalue for the operator, ap(tx), depends on the time 
derivative of the thermal Bogoliubov parameter, n+(t;p). 

Hence, we define the time dependent operators, ap{tx) and ap{tx), by 



a'^itx) =exp^-i J dtsn+^p{t,)j aP, (54) 
a"itx) = expji / ^ dt,n+^p{tA , (55) 



where we write 

nf^^it,) ^ uj+J'^P - th+{t,;p)T^^. (56) 

The time dependence of the operators, bp{tx) and bp(tx), is obtained from 
Eq. (j43|) . Time derivatives of 'qp{tx) and fj^^tx) give 

dt^Vp{t.) = -^uj.,pr^^{tx), (57) 

dt^fj;itx)^ioj^,pfl^itx). (58) 
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8 Y. Mizutani and T. Inagaki 

Applying the thermal Bogoliubov transformation, we find the time evolution equa- 
tions for bp{tx) and 6p(ta;), 

dtX^tx) = -i{iO-,p - in^{tx;p)Tor^b^p{t^), (59) 
dth;{tx) = mu)t{LU-^p - in^{U-p)To}^". (60) 

Thus the time dependent operators, bp{tx) and bp{tx), are described as 



b'^{t,)^exp\-i dUn^M\ bP, (61) 



— oo 



mt,) = bl expji y ^ dt.^^^pitA , (62) 

where the energy eigenvalue of the operators, ft°i^p{ts), is given by 

n'^J^j,{t,) = Lu^^pS^^ - in^{U;p)Tf . (63) 

Therefore the operators, and 6p, are organized into the positive and the 
negative frequency parts of the complex scalar field, 

C(x) ^ 11^^ {ap"(t.)e^P-^ + (r36p(t.)^)"e-P-^} , (64) 

C(x) ^ I {a,"(t.)e-P-^ + {bpit^rr^Te^P-^} . (65) 

Since the energy eigenvalues, VL'^p{ts), are not Hermitian, the complex scalar field 
(|64p and (j65p are not invariant under the time-reversal transformation. 
We define the canonical conjugate field. 



<(x) ^ {~^) I l^^f {«p(i.)e'P-^ - (r36p(t.)^)"e-P-^} , (66) 



^ (-Z) j |!^y^{_a^"(t,)e-P-^ + {bpit.fr.Te^P-^] , (67) 

to satisfy the equal-time canonical commutation relations, 

[C(i,a^),^f(t,y)] - [C(t,a:),vrf(i,y)] (68) 

Thus the complex scalar field is constructed by Op and bp with an ordinary canonical 
commutation relations. 

The unperturbed hat-Hamiltonian for the complex scalar field, and (^", 
should be defined to derive the equations of motion. The time derivatives of 
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Eqs. (|M )) - (p7|) give the equations of motion, 

~n+{tx;\Vx\)Tf<l>l4x)+h^{tx;\Vx\)Tfcpl_{tx), (69) 

+h+{tx; |V,|)<_(a;)ro^" - n_(t.; |V,|)<+(t.)To^", (70) 

-n+itx, iVxDTfnl^x) + h^itx, \\/x\)Tf7Tl_{x), (71) 

{dt^ + ifi)7T2{x) = -{-\7l + m')4>2{x) 

+h+{tx; |V.|)<_(x)rr - n^itx: \Vx\)nl4x)Tt, (72) 

where _(_ and 0" ^ denote the positive and the negative frequency parts of the 
complex scalar field, (IMl) and (p5|) . tt" ^ and 7f"_|_ represent its canonical conju- 
gate (|66p and (p7)) . The terms proportional to n±(t2:; IV^I) appear from the time 
dependence of the thermal Bogoliubov transformation. 

The unperturbed hat-Hamiltonian, Hq, for the complex scalar field in an equi- 
librium system consists with the kinetic term and the conserved charge as 

Ho^ J d^x[7:^{x)7:^{x) + ~^Z{x){-Vl + m2)C(x) 

-^^Ji{ra{x)K{x)-<{x)ra{x)]], (73) 

see Appendix A. In NETFD the time dependent Bogoliubov parameters induce the 
thermal counter term. 

As a simple extension of the thermal counter term for the neutral scalar field we 

set 

Q„ = / d^^':{x)^^^^^4^Tf.i{x) 

+C(x)z "-^^!' V /(-Vg. + m^)0g(x)}, (74) 

where rneans^/V^ + m? and ndtx', |Vx|) is a coefficient which consists of the 
Bogoliubov parameters. 

The time dependence of one of the Bogoliubov parameters is fixed by the self- 
consistency renormalization condition with the thermal counter term, Qn^- How- 
ever, the complex scalar field depends on two Bogoliubov parameters, ri+(i; iVxl) 
and n_(t; |Va;|). To fix all the Bogoliubov parameters, we introduce an additional 
thermal counter term. The second line in Eq. ([75)1 appears for the charged scalar 
field, i.e. no such term for the neutral field. We introduce an additional counter 
term for the second line in Eq. (|73p . 

Q-^^IJ d'x[4>:{x)fi,{tx;\\7x\)T^^7r^aitx) 

-7f,"(x)Ac(i.;|v.|)ro"^^f(x)}, (75) 
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where /ic(ix;|Vx|) is an coefficient which consists of the Bogohubov parameters. 
Thus the thermal counter term for the complex scalar field is given by 



Q — Qn Qc 



.^{C(a=)Ac(tx;|V.|)To«^f(t.) 



(76) 



The impcrturbed hat-Hamiltonian for the non-equilibrium complex scalar field, 

Hq, is represented as 

Hq = Hq -Q 

a/3 



X 



<i.x) 1 - * T Ta 7rP{x) 



/I sa0 

-7f«(x) (-^+2^0^; |V,|)Toj (a;) 



(77) 



The time dependence of the field is described by the Heisenberg equations of 
motion. Calculating the commutator between the unperturbed hat-Hamiltonian and 
the field, we obtain 



.nc{tx] \^x\) , 



'-To Trf (a;) - -Ae(tx; \V,[}T^^cp^(x), (78) 
•:ix) = Trf (x) {l - i^EM^^Toj + |V.|).^f (x)To^«,(79) 



^ (-V^ + m^)</-^(a;)--Ae(tx;|V.|)To«''7rf(a;), (80) 



{dtx +«/")< (a;) 
= -^f(x) (l- 



fe^To (-V^ + m^) + h^t.; |V.|)7ff (x)Tr. (81) 
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These equations should coincide with Eqs. (|69l) - ([7^ . This condition fixes the coef- 
ficient in the thermal counter term, 

nc(tx;\V^j:\) = , (82) 

ficit^; l^x]) = n+it^; |V:,|) - n_(i^; |V:,|). (83) 

We notice that the chemical potential, fic{tx', \'^x\) is given by the difference of the 
Bogoliubov parameters for the positive and the negative frequency parts. 



3. Self-Consistency Condition for the Complex Scalar Field 

In this section we discuss the self-consistency renormalization condition for the 
complex scalar field. For this purpose we evaluate the scalar propagator in NETFD. 
The time evolution for the field is given by the hat-Hamiltonian, H. It is divided 
into two parts, 

H^H^ + H„,u (84) 

where /fg and Hint represent the hat-Hamiltonian for a free complex scalar field 
and the interaction part in equilibrium system, respectively. 

Since the unperturbed hat-Hamiltonian for the field, (pai is not Hq but Hq in 
non-equilibrium systems, we rearrange the hat-Hamiltonian (|84|) as 

H = Hq+Hi, (85) 

where Hj is the interaction hat-Hamiltonian in NETFD, 

Hi - H^nt + Q- (86) 

The perturbative calculation should be accomplished by the unperturbed hat- 
Hamiltonian, Hq, and the interaction hat-Hamiltonian, Hj. As is shown in Ref. [1] 
the interaction hat-Hamiltonian (l86l) satisfies the condition, 



{0\Hi = 0. (87) 
In the interaction picture the full scalar propagator is given by 

D''/{t,,ty,x-y) = {0\T[^':{x)4>iiy)uioo,~^m, (88) 
where T denotes time-ordering operator and u(t,t') is the time evolution operator, 

u{t,t')=e^p (~i f dtsHiiU)] . (89) 



From Eq. (|87l) we find that the thermal vacuum, (6*1, is annihilated by the interaction 
hat-Hamiltonian and satisfies, 

{e\u{t,t') = {e\, (90) 

which is necessary to adopt the Feynman diagram procedurally. 
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Applying the thermal Bogoliubov transformation with Eq. (j90p . we rewrite the 
thermal full propagator (|88p as 

= B-\n^(^t.- IV.D) (^'ij^'^) S(rH(^,; |V,|)) 

with 

-^x)(^|<^^^_(2/)w(^y,^.)0^^+(x)^i(^„-^)|^^}, (92) 

d^,^2^^(a;,y) = 0(i, - i,)(0|0^^^.(a;)zi(t,, t,){05,+ (y)T3Pw(i„ -oo)|0) 

-ix)(^|{^«.+ (y)T3pu(t„i,)^jV(a;)7/(i,,-oo)|^), (93) 
rfc^Ca^^y) = ^(^x -i^,)(^|{r30e,-(a;)riw(t,,tj,)^^^_(y)7.(i^,-oo)|(?) 

+e{ty-u){e\^'Sy)u{ty,u){T^^^,^{x)y^u{u,-^)\e), (94) 

d^,^4^^(a;,y) = 0(t, - i,)(0|{T305,_(a;)P7.(i,, i,){^^,+ (y)T3Pu(t„ -^)|0) 
+^(ty-ix)(^|{^4,+ (y)r3pu(tj„i,){T305,_(a;)PM(i,,-w)|0}, (95) 

where ^^-i- and </>|^-|- represent the positive and the negative frequency parts of 
the complex scalar field (|44p and (j45|) . Only the diagonal elements in d]}i^ and 
d^V^ have non- vanishing values in the free propagator (jB.3| . Radiative corrections 
generally induce all the elements in Eqs. (|92|) - (|95l) . 

The time evolution operator u{t,t') is expanded perturbatively in terms of the 
thermal counter term, Q and the interaction part Hint- First we evaluate the con- 
tribution in the leading order of the thermal counter term, Q, at the tree level with 
respect to Hint- Performing the thermal Bogoliubov transformations in Eq. (|76p . 
we rewrite the thermal counter term, Q, as 

Q = ~'J {^'Hitx;p)il(l + n^{t.;p)vWp} ■ (96) 



Substituting Eq. (|89|) and Eq. ([96|) into Eq. (|9T|) . we obtain in the leading order of 
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the thermal counter term, Q, 



(Pp 1 
(27r)3 2^^' 

d{tx-ty) I dtsfi+its^p) + 9{ty ~ tx) I dtsh+{ts;p) 



X B-\n+{t,;p)) (^|jJ^B(rH(i,;p)) 

J^P_J_^iuj^.j,{t^-ty)^-ip-{x-y) 
(27r)3 2wp 



X {6{t^-ty) / dt^n_{ts]p)+d{ty-t^) I dtsn_(t^;p)\ 

\ J — oo J —oo j 

X T^B{n^{tx-p)f (^J B-\nAty-p)fT^. (97) 



These counter terms contribute d}^^ and 4 in Eq. ([9T|). 

Next we evaluate the contribution from the interaction, llint- The contribution 
to the perturbed propagator is represented as 

y" d^zid'^z^D^^' , X - (^^^ ^ ^ _ Z2)D''^'\u, ,ty,Z2-y) 

= Jdt^.dt,, [B-\n+{tx;\Vx\)r^'Si:l'J'{x,zi,Z2,y)B{n+{ty;\^y\)r-^ 

+i?-i(n+(t,;|V.|))"^ME^^J(x,zi,Z2,2/){S-^(n-(t,;|V,|)fr3p^ 
+{r3B(n_(t,;|V,|))^}"-'i5I];^^J(a;,zi,Z2,y)i?(n+(t,;|V^|)r^ 
+{r3B(n_(i,;lV,|))^}"^^JI];^^J(a;,zi,Z2,y){i?-^(n-(t,;|V,|))^T3p^ (98) 

where D^^ represents the free propagator (jB.3p and 

^E^^J(x,zi,z2,y) (99) 

^ f D^ljiix - zi)i^R{t^^,t^^, zi - Z2)Dli1jf {z2 - y) 5T.]^-^^{x,zi,z2,y)\ 
{ Dl%{x-zi)iJ:A{U.„t,„zi-Z2)Dl%{z2~y)J ' 

<5I]I,^J(x,zi,z2,y) (100) 

'5E^i2(2;'^i'^2,y) - DI^r{x - zi)iT.R{t^^,t^^,Zi - 22)1)^1^^(22 - yY 



'Dl^ix ~ zi)iJ:A{tz,,U2,zi -Z2)Dl]^{z2 -y) 
<5E^^J(x,zi,Z2,y) (101) 
- Dl\{x - zi)iY.Aitzi,tz2:Zi - Z2)Dl^p.{z2 - yY 

-Df^Ai^ - ziY^R{tz,,tz2,Zi - Z2)Dl]j^{z2 - y) (5i;|'^3(x,zi,Z2,y) 
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5Y?^J{x,z^,Z2.y) (102) 
\5T?j^^^{x,zx,z2,y) Dl\{x - zi)iT.R{tz^,tz2,zi - Z2)Df^^{z2 - y) ^ 



with 

^12 



Z2 



i(a;,zi,Z2,y) 

= Dl]j^{x - zi){iY}'^{tz^,tz2,Zi - Z2) + iY.R{tzi,tz2,Zi - Z2)n+{tz^; jV^ 

-n+{tz, ; I V,, , , z^ - Z2)}Dl%iz2 - y), (103) 

^5^5,2(2;, 2^1, 2^2, y) 

= DI]r{x - zi){ij:^^{tz,,tz2,zi - Z2) +i^R{tz,,tz:,,Zi - Z2)n-{tz2; iV^J) 
+n+{tz,;\\/zA)i^A{tz,,U,,Zi - Z2)}dIIa{z2 - y), (104) 

^^5,3(2^' ^2, y) 

= --C)g!4(a^ - zi){iT,'^'^{tz,,tz2,Zi - Z2) + i^R{tz,,tz2,Zi - 2;2)n+(t^2; 1^22!) 

+n_ (t,, ; I V,, |)zSA(t.i , , - ^2) }i?o'fl(^2 - y), (105) 
'5^3,4(0;, zi,Z2,y) 

= D^;^^(2; - zi){-iS2i(i,^ , tz,,zi - Z2) - i'SRitz, , tz, , 2:1 - Z2)n^{tz, ; | V;,, |) 
+n^{tz,;\Vz,\)t^A{tz„tz„Zi - Z2)}Dl\{z2 - y). (106) 

The self-energy, S"*^ , has a 2 x 2 matrix form in the thermal doublet notation. The 
retarded and the advanced parts of the self-energy, E/j and T^a, are defined by 



The first and the last terms in the right-hand side of Eq. (|98)) have the same Bo- 
goliubov transformation structure as the first and the last terms in the right-hand 
side of Eq. ([W)) . respectively. 

In Ref . 121] the self-consistency renormalization condition has been proposed for 
a non-relativistic field by imposing {O\^fuU,p{tx)^fuU,k{tx)\0) = 0, where the sub- 
script "full" denotes the perturbed operator with all order of radiative corrections. 
In Ref. 1^ it has been shown that the self-consistency renormalization condition 
induces the quantum Boltzmann equation for a neutral relativistic scalar field. The 
self-consistency renormalization condition is naturally generalized for the complex 
scalar field, {9\^fuii,p{tx)ifuii,k{tx)\6) = and (0\fifuH,p{tx)VfuU.k{ix)\0) = 0. The 
conditions fix the terms, dl\ and in Eq. (PT|) H Thus we obtain integral equations 

^ Substituting the fields Il44|l - lj45|l and taking the equal time limit, we obtain 

lim dll(x,y) = - [ f ^3 j_ ^ e'P-^e-'fe-3^(e|T[g,(f.)4ft.)»(c^,^^)P), 

lim dll,{x,y) = - f f P f ^ L_ ^ e-^P-''e^''-y{e\T[ij,{tMt.Moo,-^)]\e). 
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for n±(t;p), 



(108) 

/ dts [ ^-^e-'-P '-'^-y^ ^h^{ts;p) + YiuY [ dt,,dt,JJ:l\{x, zi, Z2,y) ^ 0. 

(109) 

After the spatial Fourier transformation the self-consistency renormalization condi- 
tions are rewritten in the t-representation. Acting the time derivative operator, dt^ , 
on Eqs. (fTOS]) and (fTog]) . we find 



(110) 



h+{tx\p) = -2wp5t^|^lim j dtz^dt^JY}^^^{t^,,tz^,tz^,ty]p)^, 
n_(ia:;p) = -2wp(9t, l^hm j dt:,^dt:,^6Y.'^B^{tx,t^_^,t^_^,ty;p)^, (111) 
where 5Yjb,i and 5YiB,i are given by 

= DaMiix -t,,]p)^iT,^^(t,^,t,^;p) +n+{t,^;p) iT,R{t,,,t,^;p) 

~n+{t,,-p) iJ:A{t,,,t,,;p)}Dl%{t,, ~ty;p), (112) 

= Df,Ait^ -t^^;p)^~iY.'^^{t^^,t^^;p) - n^{t^^;p) iY.R{t^^,t^^;p) 

+n-{tz,;p) iT,A{tzi,tz2;p)]DllA{t^^ - ty;p). (113) 

The equations (jllOp and (jllip describe the time evolution of the thermal Bogoli- 
ubov parameters, n±{tx;p)- In a similar manner developed in Ref. [22] it can be 
shown that the Bogoliubov parameters, n±{txTp), for complex scalar field cor- 
respond to the particle and anti-particle number densities under the self-consistency 
renormalization conditions ([T08| and (fTOQ]).^ ^'^ Thus Eqs. (fTTOl) and (fTTT|) give 
the quantum Boltzmann equations for the charged Bosons. 



4. Boltzmann Equation for the Complex Scalar Field 

The time evolution equations for the Bogoliubov parameters (jllOp and (jllip rep- 
resent the transport equations for the number densities. In this section we discuss 
the validity of the time evolution equation. For this purpose we perturbatively eval- 
uate Eqs. (jllOp and (jllip for interacting complex scalar field. One of the simplest 
interactions is a four-point interaction between complex scalars. We start from the 
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Hamiltonian with a A((/)jj(/>a)^ interaction, 



{nl{x)na{x) + 0l(x)(-V2 + m'')<Pa{x)] 



'i^l {(l)\{x)T:a{x) - ■Kl{x)(j)a{x)] + -(l)\{xf (j)a{xf 



(114) 



In NETFD the Hamihonian (|114p is extended to the hat-Hamiltonian which is 
defined by 



H 



d^x 



{~4>Z{x)K{x) - K{xWa{x)] + 5] - {{Ux)r,V^l{x)Y . (115) 



7=1 



Employing the Feynman rules, we calculate the off-diagonal elements (|103p and 
()106|) in the self-energy. In the thermal doublet notation the propagator (jB.SP is 
assigned for an internal line as Dq^t^^ . For a vertex we assign, (— i)A", which is 
given by 



(-z)A" ^ (-z)A 



-1 



(116) 



Since the radiative correction for the propagator has no momentum transfer to 
the internal line at the one-loop level, the self-energy has a diagonal form up to 
the one-loop level. The time evolution of the thermal Bogoliubov parameter can be 
induced from the 2-loop self-energy illustrated in Fig. [TJ We compute the diagram 
and find 

x{Do{t,, , t,, ; ki)T3r^{Doit,, , t,, ; fc2)T3r^(-i)A^{i?o(t.. , i.^; fe3)T3}''^ (117) 




Fig. 1. 2-loop thermal self-energy in \{<f>a)-^ (j>^ interaction model 
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After the contraction of the index "7" Eq. (|117p readt0 

11=112=113=1"' V ; ^ ; ^ ; 



)3 SWfeiWfcjWfeg 

x(27r)''^(5(3'(p - fci - fc2 + fc3)e-*(^^''=i''i+^<=^''=2.-2--Bc,.3..3)(*-i-*-2) 

Fiitz2;ki, ^2, ^3, «1, «2, «3) -F2{tz2;ki,k2, h, il, «2, «3) 



Fi{tz^;ki, k2, ks, 11,12,13) --^2(^22;^!) ^2, ^3, «i, «2, i3) 
F2{tzt;ki, k2,k3,ii,i2, is) --^2(^21; fci, ^2, fcs, ii, 12, 13) 



,(118) 



,-Fi(izi; A:i, k2,k3,ii,i2, is) -Fi{tz^;ki,k2, ks,ii, «2, 13) , 
where we write as 

Fi{t;ki,k2,k3,ii,i2,is) = Iii{t;ki)fi^{t;k2)Jt:,{t;ks), (119) 
F2{t;ki,k2,k3,ii,i2,is) = fi^it; ki)fl^{t; k2)fi3{t; ks), (120) 

and 

Ec,q,l = ^^+,q = '^g - Ec^q^2 = -^^-,9 = -^^g ^ (121) 

/i(t;q) = n+(i;g), /2(t; g) = 1 + n_(t; g), (122) 
h{t;q) = l + n+{t;q), f2{t; q) = n^{t; q). (123) 

We substitute the self-energy (jllSp into Eqs. (jllOp and (jllip and integrate over 
the time variables. Taking the equal time limit, we obtain 



lim / dt^-^dt^^6'E^B,i{tx,tzi,tz2,ty;p) 



*^ 1 A2 ^ ^ ^ /■ rf3fei d^k2 d^k3 1 
_^ ^^^^Y f^J^ J (27r)3 (27r)3 (27r)3 8c.pC.,iC.fe2^fc3 

Sin{(w+ p - E'er ji - fe2,i2 + -E'c,fc3,i3)(^a; - is)} ,„ ^3^(3), , , , , n 

X ^ ^ — '-{2T:)-'S'--'>ip-ki-k2 + k3) 

t^+,p - J^C.kiAi - J^C,k2.i2 + ^C.ka.ia 

x{n+{ts;p)U^{ts;ki)fi^ (ts; fc2)/i3(is; ^3) 

-(1 + n+{t,;p))f,, iU; ki)f,, {U; fc2)4(t,; ks)} , (124) 



lim / dt^^dty_^6Yi^B^{tx,t^_^,tz^,ty;p) 



1 /■ d^ki d^k2 d^ks 



°o 2 ^fr'i.fr'i.fr'iy (27r)3 (27r)3 (27r)3 Sup^.^uJk^uJk, 

^ sin{(c.„,, + EcM.. + EcM,. - Ec,k,,i,){t. - U)] (^^)3^(3)(p _ fc^ _ fc, + ^3) 

UJ-^p + Ec,ki,li + Ec,k2,l2 ~ Ec,k3,l3 

x{n_{ts;p)fi^{ts; ki)},^ {1^; fc2)/j3(is; ^3) 

-{l + n^{U;p))hAts-M)k{ts-M)h,{ts;k3)]- (125) 



''A simpler formalism is developed to calculate the self-energy in Ref s . [28l and l29l 
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Though the chemical potential seems to contribute the frequency and the amplitude 
of the time oscillation, Eq. (|125p is independent on /i because of the relationship 
(|T2T]1 . 

Thus we find the time evolution equations (|110l) and (|llll) for the thermal Bo- 
goliubov parameters of the complex scalar field, 

A2 , f d^ki d^k2 d^ks 1 



A - f 

n+{t^;p) = -— XI XI XI / 



il =1 '2 = 1 *3 = 1 ' 



(27r)3 (27r)3 (27r)3 SwpWfei Wfc^ 



xcos{(wp - Ek,,^, - Ek.,,^^ + Ek^^i^){t^ - ts)} {2n)^ 5^^'> {p - ki - fca + fcg) 

-{l + n+{ts;p))hAts:ki)U,{ts-M)k{ts-M)], (126) 
A2 ^ ^ ^ „ [ d^ki d'^k^ d^ks 1 



(27r)3 (27r)3 (27r)3 SwpWfciWfe^Wfeg 



'^-(^-P) = -yE EE '^^^ /■ 

^l=1^2=li3=l•^"°° 

XCOs{(Wp + Ek,A^ + Ek.,.i^ - Ek,,^,)it^ - t s)}i2Tr fS'^^^p - fci - fcs + feg) 
X {n_{ts;p)f,, {ts; ki)f^^ {1,,; k2)fi3 [ts] fcs) 

^{l + n-{U-p))f,,{UM)Mts-M)hAts;k:,)], (127) 

where we use 

Eqs = ^g, Eq^2 = -C^g. (128) 

These equations describe the thermal evolution of the Bogoliubov parameters for the 
complex scalar field with the \{(j)\(l)aY interaction. It has the consistent statistical 
structure with the quantum Boltzmann equations for the two body charged Bose 
particle scattering process. 



5. Canonical Quantization for a Dirac Field 

We quantize the Dirac field in a similar manner to that used for the complex scalar 
field developed in the previous sections. Because it is easy to distinguish, we use 
the same expressions for the Dirac field with those for the scalar field, for example, 
the creation and annihilation operators, mass and so on. 

We review the framework of NETFD for the Dirac field, simply. Ordinary and 
tilde operators for the Dirac field obey the anticommutation rules, 

{a^,a^J} = {a^al^} = {2nf5^^\p - k)S^-\ (129) 

{b;X^} = H^K^} - (2^)^<5(3)(p _ fc)5-, (130) 

where the subscripts r and s represent two-component spinors. All other anticom- 
mutators vanish. We use the tilde conjugation rules given in Eqs. (I4])-(I7]). For the 
Dirac field it is convenient to use the thermal doublet notation defined by 

"^" = (z|t)' s;- = (a;t-,a;). (131) 
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It should be noted that the definition is different from Eq. ([T| 

Time-dependent thermal Bogoliubov transformations are introduced for each 
operator with two independent thermal Bogoliubov parameters, n±(t;p)^ 

i;'^{t) = B{n+{t;p)r"af{t), Cp^it) = af {t)B-\n+{t;p)f^ , (132) 

r^l'^t) = B{n^{t-p)r^bf{t), n^it) = bf{t)B-\n^{t-p)f'^, (133) 

where the thermal Bogoliubov matrices are given by 

(134) 

B-\n^) = (\ -''t V (135) 



1 1 - n± ^ 

Here the thermal Bogoliubov parameters are assumed to be independent on the 
spinor index, s. Thus the thermal Bogoliubov parameters depend only on time, t, 
and the absolute value of the momentum, p, in the homogeneous and out of equi- 
librium system. The transformed operators, and 77^, and their tilde conjugates, 
and r^p, obey the anticommutation relations, 

= {tr>,£,l^} = (2^)3<5(3)(p- fc)5-, (136) 

H^r^l^} = = {2irf5^^\p~k)5^^, (137) 

with all other anticommutators vanish. The thermal vacuum, \6), is defined to be 
the state such that 

^pl^)=l;i^)='7;i^)='?;i^)=0, (138) 

mi} ^ {e\t^ = {OW^ ^ {ew^ = G. (139) 

The thermal expectation value is obtained by the expectation value under the ther- 
mal vacuum, \0). 

Below we quantize the Dirac field in the interaction picture. To keep the time 
independence of the thermal vacuum, the Dirac field has to be decomposed as a 
Hermitian form with the transformed operators, 

Ux)^ J ^^^{C^(i,)«^(p)e'P-^+r?;t(,^),.(p),-.p-a=}, (140) 



with 



and 



C;{t.) = e;e-^"+-*% CpHt,) = (142) 

v;it.) = ?7;e— rj;\t,) = 77^te--..*==, (143) 

w+,p =ujp- (144) 

uj^^p =ujp + 11, (145) 
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where uip is the relativistic energy eigenvalue. The chemical potential, /i, is intro- 
duced according to the discussion in Appendix A. The functions, M-'(p) and v''{p), 
are the solutions of the Dirac equation with the positive and the negative frequency 
parts, respectively. We write = u^^'y^ and = v^'^j^. The eigenfunctions obey 

(7V-7-p-mK(p) =0, (146) 

(-7°^^ + 7 • P - m)v' (p) = 0, (147) 

u^(p)(7%p-7-p-m) =0, (148) 

v'{p){--f°ujp + 'r -p-m) = 0, (149) 

and 



7°Wp - 7 • p + TO, 



S 

s 

The tilde conjugate for the Dirac field is given by 
d^p 1 



TO. 



(150) 
(151) 



^^^'^^^ / (t^3 E {m^X'(P)e-'P-'' + %\t,y*{p)e^P-''} , (152) 



r)3 ^/2il^ 

(2^771^ 



The time dependence of the non-tilde and tilde operators is combined into the 
thermal doublet notation, 



The Dirac field is rewritten in the thermal doublet form, 
d^p 1 



(154) 
(155) 



(27r)3 



J (27r 



d^p 1 



d^p 1 



E {ir(i-)^'iPy^''' + {r3%{t.frv^{p)e-^P''} , (156) 



+ 



(27r)3 

Uit.) -ivp{t.))v%p)e'P-^} 
d^p 1 



= I ^-i= 5: {irit.)u^{p)e-'P-^ + {v;{t.frsrv%p)e'P-'^} • (157) 



These fields satisfy the equal-time anticommutation relations, 
{4''{t,x),i^l{t,y)} = S^'\x-y)6'^^, 
where we write Vj" = V'|^7°. 



(158) 
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Next we decompose the Dirac field in terms of the original operators, a* and 
bp. The time dependence of the original operators is fixed by the time evolution 
equations for the transformed operators. The time derivatives of Eq. (jl54p give 

dtj'p"it.) - -iLo+.p^p^it^), (159) 
dtj;-"it,)^iuj+j''''{t,). (160) 

Applying the thermal Bogoliubov transformation (jl32|) with Eqs. (jl59p and (|160p . 
we obtain the time evolution equations for the original operators, ap"(ta,) and 

dtXp^it.) = -i{uj+.p + ih+{t^;p)Tor^af{t,), (161) 
dt^K'^i^-) = afit,)i{^+,p + in+{U;p)Tof", (162) 

where Tl^^ is given by Eq. ([5^. Eqs. (|16ip and (|162p indicate that the energy 
eigenvalue for the operators, ap°'{tx) and ap°'{tj;), depends on the time derivative 
of the thermal Bogoliubov parameter, njf.{t;p). The solutions for these equations 
are found to be 

a;'"(i,) =exp|-iy dtsn+,p{ts)[ af, (163) 
ar(ir.) = af exp |z J ^ dt,r!+,p(i,)| , (164) 

with 

f)fp(i,) = i.+^pS'^^ + ih+{ts;p)Tf . (165) 

The time dependence of the operators, bp'°'{tx) and bp'°'{tx), is found from 
Eq. (|155[) . The equations of motion for rjp^itx) and fjp'^itx) are given by 

dt^v'^^'itx) = -zc^-.pry^'"(<.), (166) 
9t^v''''{tx) = ic^-,pf;^'"(i.). (167) 

Applying the thermal Bogoliubov transformation, we find the time evolution equa- 
tions of the operators, bp'°'{tx) and bp°'{tx), 

dtj;'"{tx) = ~i{^~..p + ih^{tx;p)Tor^bf{tx), (168) 

dtip"{tx) = bf{tx)i{oj-,p + in-{tx;p)To}''". (169) 

The solutions for these equations are represented as 

b;'"{tx) = exp j-^y ' dt,f]_,p(t,)| bf, (170) 

b;'"{tx) = bfexp |i I ' dtsn^M^ , (171) 

where 

nl^pits) = u;-,pS"^ + ih^{ts;p)Tf . (172) 
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Hence the Dirac field is decomposed in terms of the original operators 
(fp 1 



ra{x) 



(27r)3 yf^Tp 

dPp 1 
(27r)3 J2^p 



S 

(173) 

S 

(174) 



It is easy to find that tliese fields satisfy the equal-time anticommutation relations, 

{V4"(i, x), (i, y)} = S^^Hx - y)5'^P, (175) 

where we define i/'i" = i^al^ ■ 

Time evolution of the Dirac field is generated by the hat-Hamiltonian (|13p . The 
unperturbed hat-Hamiltonian, iJg, for the Dirac field, "0" and V'" j t)e found 
from the equations of motion for the Dirac field. Differentiating Eqs. (|173p and (|174p 
with respect to time and using Eqs. (|146p - (|149p . we derive the equations of motion 
for the Dirac field, 

^I'dtM^) = Hi • V, + rn - M7°)C(a^) 

+z7V(i.; |V.|)To"^v£+(a:) - i^'>h^{t,; |V,|)To"'5<_(x), (176) 

-^dtj^{xh° = C(a;)(»7 • V, + m - fi^") 

|V.|)V;f,_(a:)7°rr " *n_(t.; \S7,\)i;l^{xh°Tt , (177) 

where V'a,± ^^.d V'a.i describe the positive and the negative frequency parts of the 
Dirac field (|173p and (|174p . respectively. 

The unperturbed hat-Hamiltonian in an equilibrium system for the Dirac field 
at finite chemical potential is given by 



(178) 



We introduce the thermal counter term for the Dirac field in the following form. 



--C(^)/i.(i.;lv.|)7X^^f(x) 



(179) 



where nd{tx', \'^x\) and fid{tx] \'^x\) are the coefficients consist with the Bogoliubov 
parameters, n+{tx] \'^x\) and n-{tx] |Vx|). The unperturbed hat-Hamiltonian for 
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the non-equilibrium Dirac field, iJg, is defined by 
Hq = Ho - Q 



C(^) ( 1 + ^^iM^ToV' (-^7 • V. + m)^f (x) 



(180) 



The Heisenberg equations should reproduce the equations of motion (|176p and (I177P . 
This condition fixes two types of the coefficients, nd{tx] |Va;|) and ^d{tx] |Va;|), in 
the thermal counter term as 

(id{tx\ |V,|) = h+{tx] |V,|) - |V,|). (182) 



These equations have the same form with those for the charged scalar field ((82|) and 
(El. 



6. Self-Consistency Condition for the Dirac Field 

In this section we introduce a self-consistency renormalization condition for the 
Dirac field with an interaction. Hint- As is shown in the previous section, the un- 
perturbed hat-Hamiltonian (|180p contains the thermal counter term in NETFD. 
In a similar manner with §3, the interaction hat-Hamiltonian, Hi, is given by 
Hint + Q- The perturbative calculation should be developed under the unperturbed 
hat-Hamiltonian, Hq, with the interaction hat-Hamiltonian, Hi. 

6.1. The Bogoliubov Parameter and the Number Density 

Following the procedure in Ref. [22l we show that the Bogoliubov parameters for 
the Dirac field correspond to the number densities under the self-consistency renor- 
malization conditions. In TFD the observable particle number densities are given 
by the vacuum expectation values of the number operators, 

(2^)^5-5(3) (p - k)nHAt;p) = {9\a%iLi,muiLkim, (183) 

(2^)35-5(3)(p _ k)nHAt:p) = mb^Upm^ii^i^im, m) 

where the lower index "/wZ/" indicates the perturbed operators. 

We define the Bogoliubov transformation for the perturbed operators, 

e/:ujt) = 5("+(i;p))"^a}f,,p(i), ^}ZiJt) = -aX,,Jt)B-\n+{t;p)r'^, (185) 

VfZiJt) = S(n-(i;p))"^<„,,(i), fj}':uA*^ = 6}f,,^(t)B-i(n_(t;p))^", (186) 

with the Bogoliubov parameters, 'n±{t;p). From the condition, {0\Hi = 0, it is found 
that the transformed operators, ^f^uit) and satisfy 

mfuu,pit) = mfuuA^) = 0, efun,pm) + o, lu^^m) + o. (18?) 
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Substituting Eqs. and (IW|) into Eq. (IT55)) . we obtain 

= mtfunAt)^UiMo) + n+{t- kmtf^uA^yCfiii^Me). (iss) 

This equation implies 

(2^)35-5(3)(p _ - n+{t-p)] = (189) 

Thus the Bogohubov parameter, n^{t;p), coincides with the number density, 

under the condition, {0\^fuu,p(t)CfuU,k(t)\^) = 0- 
We also derive the condition for the antiparticle number density (|184p . The 
operators, VfuH,pit) and Vflu.pit)^ satisfy 

{(^iVfUpit) = m}Upit) = 0, v}uu,pm) + 0' %uii,pm) + «• (i^o) 

In a similar manner to Eq. (|188l) we obtain 

(2^)35-5(3) (p _ fc) {nn,- (t; V) - n- {t; p)} ^ t (t) 1^) • (191) 

Therefore the Bogoliubov parameter, n-(t;p), also coincides with the anti-particle 
number density, nH-{t;p), under the condition, {d\flfuiip{t)'>lfuiik('^)\^) = 0- Thus 
we impose the self-consistency renormalization conditions for the Dirac field, 

{muii,pmu,km) = (192) 

6.2. Time Evolution for the Bogoliubov Parameters 

The time evolution is derived from the self-consistency renormalization conditions. 
Here we begin with the perturbative expansion for the full Dirac propagator in 
NETFD. In the interaction picture the full Dirac propagator is written as 

S^^{t,,ty,x~y) = ((?|r[C(x)V5f(2/)^(oo,-^)]|(?), (194) 

where the operator u{t,t') has the same form with Eq. Applying the thermal 
Bogoliubov transformation, the full propagator (|194p is rewritten as 

(.-(L) f (x!^)) ^("^^^-'^^'^^ 

+T,B{n^{t,-AWx\)r [^'^^^y^^ ^J^^^^^ (195) 
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with 

sT''^x,y) = e{t,-ty){e\4,l^{x)u{t,,ty)4,l-_{y)u{ty,-^)\e) 

-t.)(e|V5^,^_(2;)?^(t„i.)^J,V(a;)ii(i.,-c^)|^), (196) 
sY'-{x,y) = -i,)(0|V'5^_V(a:)^.(t,,i,){V;?,+(2;)T3Pu(i„-oo)|0) 

s2^^^(a;,y) = -i„)(0|{T3V'5,~(a:)pu(t.,t,)^J^^(2;)M(t^,-(X))|0) 

-t,)(0|V;j.^(2;)M(ty,t,){r3V'j,_(x)pu(i„-w)|0), (198) 
sp^=(a;,y)=0(t,-i,)(0|{r3^4,-(a;)r^w(t,,t^){V:;e,+ (y)r3pw(t^,-^)|0) 
-0(t^-t,)(0|{V;4,+ (y)r3P?/(iy,t,){T3^5,_(x)p?/(i,,-oo)|0), (199) 

where ■0|'± 8'iid ^l^-i- show the positive and the negative frequency parts of the Dirac 
field, (|156p and (|157|) . respectively. 

We expand the time evolution operator u{t,t') in terms of the interaction hat- 
Hamiltonian (I86p and evaluate the propagator in the leading order of the thermal 
counter term, Q. After the Bogoliubov transformation the thermal counter term, 
Q, is written in terms of the transformed operators, and 77^, 

Q = - j (0jE{"+(^-^')^p^^p^-'^-(*-^^)^pV}- (200) 

Substituting Eq. pOOp with Eq. (l89l) into Eq. (I195P and dropping the higher order 
corrections, we obtain 

{e\T[i:2{x){pi{y)u{^,-^)]\e) - {e\T[ra{xWa{y)]\e) 

(27r)3 2ujp 



xi?-i(n+(t,;p))(^[jjji?(n+(t,;p)) 
(27r)3 2wp 



X \9{tx-ty) I dtsn^{ts;p) + 0{ty -U) / dtsn_(is;p) 

y J —00 ^ —00 / 

X T3B{n^{t^;p)f (^J B-\n^{ty;p)fT3. (201) 

A finite correction appears for s}^ and in Eq. (jl95p . 

Next we consider the radiative correction from the interaction. Hint- At the 
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leading order it is represented as 

= yd4zid4z2[i?-i(n+(t,;|V,|))"^MI]];;f (x,zi,Z2,2/)S(n+(t,;|V^|))^^^ 

+B-\n+{t,;\\/,\)r'''SJ:l\J'{x,zi,Z2,y){B-\n^{ty;\^y\)fTsr^^ 
+{r3i?(n_(t,;|V,|)m^5E],;3^^(x,zi,Z2,2/)S(n+(t,;|V,|)r'5 
+{T3i?(n_(t,;|V,|)rn^5I]J,;2'(:.,zi,Z2,y){i3-i(n_(t,;|V,|)rr3p^], (202) 
with 

(5E^;f (a;,zi,z2,y) (203) 
^ [ Soljfix - zi){-i)^R{t^^,t^^, zi - Z2)S^lj^{tz^ -ty) 5i:f-^^{x,zi,z2,y)\ 

V ' Sl%{x - zi){~i)^A{U„U,,zi - Z2)Sl%{z2 - y) y ' 
5^l';'{x,z^,z2,y) (204) 

^ / 5Y}j^2{^,zi,z2,y) S^Ij^{x - zi)iJ:ji{t^^,tz^, zi - Z2)S^'^a{z2 - y)' 
~ V'5'o,^i?(2;-2i)iSA(i^i,t22,^i -2;2)5'o"a(2:2 -y) 
5T,]^J^ {x,zi,z2,y) 

S'o"^(a; - zi)iEA(i2i,iz2,2i - 2;2)5'^,^i^(22 - y)' 

\'5'o,a(2; - 2:i)«S/?(t^i,t^2,^i " ^2)S'o,i?(22 -y) (5i;|?3(a;, 2;i,z2,y) 
5E;^^2'(a^:^i.^2,y) (206) 

^fSl^y,{x-zi){~i)^A{t.,,t,,,zi-Z2)Sl]A{z2-y) \ 

V 5T.f,i{x, zi,z2, y) 5*^ a(x - zi){~i)T.A{t:,^ , i^^, - 22)5^,^(^2 - y) / ' 

where S'q is the free propagator given by (IC.3P and 

5Y}p^-^{x,zi,z2,y) 

= SI]r{x - zi){~iY}'^{t^^,t^,,_,Zi - Z2) " {-i)T.R{t^^,t^^,Zi - Z2)n+{t^^;\^ ^^\) 
+n+{t,,;\V,,\){~i)^A{U,.U,,zi-Z2)]Sl%{z2-y), (207) 
(5E^yx, zi,z2,y) 

= SI]r{x - zi){-il]"(t2i,i^2,2;i - Z2) +iEfl,(izi,iz2J^i - ^2)?^-(^^2; I^^J) 

+n+{t,, ; I V,J)iS^(i,, , i,,, 21 - Z2)}SI]a{z2 - y), (208) 
<5SF^3(a;,zi,Z2,y) 

= '5'oiA(a; - 2;i){jS22(t^i,t22,^i - Z2) + {-i)J:R{tz^,t^^,Zi - Z2)n+{t^^; \^ zj) 

+n_ {t,, ; I V, J ) (-z)S^ {tz, , t,, , 21 - 22) } (^2 - y ) , (209) 
(5E^^4(a;, zi,Z2,y) 

= 50,^(2; - 2;i){iE2^(t2,,i^2,2i - Z2) + {-i)T.Ritz,,tz^,Zi - 22)"-- (i^2 ; |^^2 1) 
-n-(i.i ; I V,J)(-z)S^(t,, , tz, , 21 - 2;2)}4i^(z2 - y). (210) 
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The retarded and the advanced parts, S/f and E^, are defined by 

Sfl, = + = + 2^2, = - = - (211) 

The self-consistency renormahzation conditions, (6'|C/t,;/,p(*a;)C/«;/,fc(*2:)|^) = 
and {6\''ifuii,p{'tx)v}uU.ki^x)\0) = 0, can be satisfied by fixing the terms, sj^ and 
in Eq. (flM)) at the equal time limit@ From Eqs. ((20T|) and ((202)) we can find the 
self-consistency renormahzation conditions for the Dirac field, 

lim d'^zid'^Z2ST,]?^^{x,zi,Z2,y) ^0, (212) 



^ ' J (2^)3 2a;p ^ ^-^^s^P) 

lim [ d^zid^Z26J:f^^ix,zi,Z2,y) =0. (213) 

The t-representation is more convenient in practical calculations. Performing the 
spatial Fourier transformation and acting the time differential operator, dt^ , on 
Eqs. (|2T2|) and ([2T3l) . we obtain 

^ h+{tx;p) = dt^ lim I dt^^dt^^5T}p^^{tx,t^^,tz^,ty]p) = 

(214) 

^ ?i_(ij;;p) = -5f lini j dt^^dt,_JT.'^p^^{t,j,,t^^,t,_.^,ty]p) ^ 

(215) 

with 

+n+{t,,;p){-i)^A{t^,M^;p)]Sl%{U,-ty;p), (216) 

Substituting the fields II156I I-I I157I I and taking the equal time limit, we obtain 

lim sf{x,y) = -^^ ^ ^ ^ ^ ^ ^ (p)n'-(fc)e>P-e-'=-2/ 

x(e|rK^(t,)c^(i,)«{c»,-^)]|e>, 

lim sf{x,y) = -^f ^ ^ ^ ^ ^ ^ (p^ (fc)e-P-e-'=-3^ 



(27r)3 (27r)3y2i;;}^^2Zj^ 

x(e|TR;(i,)r,;;(t,)«(oo,-oo)]|e> 



s r 
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-n^{tz,-p){~i)^A{tz,,tz,-p)]Sl]A{tz,~ty;p). (217) 

As is shown in §6.1, the BogoHubov parameters, n±{tx',p), under the self-consistency 
renormahzation conditions (I214p and (j215p correspond to the particle and antipar- 
ticle number densities. Thus we conclude that the time evolution for the particle 
and anti-particle number densities is described by Eqs. <\21A\ and (|215l) . 



7. Boltzmann Equation for the Scalar and the Dirac Fields 

Here we consider a Yukawa type interaction between neutral scalar and Dirac fields. 
To show the validity of the time evolution equations (I214p and (|215p . The time 
evolution of the thermal Bogoliubov parameters is studied for the scalar and the 
Dirac fields. We start from the Hamiltonian, 



H = (Px 



ipaix){-ij ■ + nip - HJ°)^a{x) 



+ ^{Traixy + Mx)i-'^l+m%)(t)a{x)}+9'4'a{x)Mx)Mx) , (218) 



where 4>a describe the neutral scalar field and g is the Yukawa coupling constant. In 
this section all the quantities with subscripts B and F represent ones for the scalar 
and the Dirac fields, respectively. In NETFD the Hamiltonian (|218p is extended to 
the hat-Hamiltonian ([T^ . 



H 



d^x 



C(a;)(-»7-V.-Hmi.-/i7")CW 



+ Y.9''{i^a{x)T^Vra{^)4>l{x) 
7=1 



(219) 



where the neutral scalar field in the thermal doublet notation, and c/)^ , and their 
canonical conjugate, tt^ and tt", are given in Ref. [22l 

In Feynman diagrams the thermal propagator, Sq^t^'^ , is assigned to the inter- 
nal fermion lines. (Se e A ppendix C.) The free thermal propagator for the neutral 
scalar field is given bjE21 



Dfit,,ty;p) = B-'{nBit,;p)r'D^,^^^it,-ty;p)B{nB{ty;p)r^^ 
+{r3BinBitx;p)fr'''D^,]2"ito.-ty;p){B-'inBity;p)fT,y^^, (220) 
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where we write 

s 1 





-ty;p) 




-ty;p) 








-y) = - 



2ujp 



ZLUp (^t^ ty ^ 

1 



2a;p ' 

/ ^^_„^"p(tx-tH) 

^^2c./ 



(221) 
(222) 
(223) 
(224) 



other components — 0. 



We assign D^^t^^ on the internal scalar lines. The Feynman rule for the vertex is 
given by, 

H)5^ = H)5(_\), (225) 



where 7 denotes the thermal index. 

We evaluate the time evolution equation for the Bogoliubov parameters at the 
1-loop level. The self-energy for the neutral scalar field is illustrated in Fig.[2][a). It 
is calculated as 

X (-l)r3"^tr [{-i)g^So{t., , ; fci)T3r^(-z)/{5o(i., , ; k^Y^f-^] ,(226) 
where tr[- • • ] shows trace manipulation on the spinor indices. Contracting the ther- 




Fig. 2. l-loop thermal self-energy derived three point scalar interaction model 
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mal index "7", we obtain 

2 2 

il = l 22 = 1 

EF,ki,ii^F,k2,i2 — ^1 • ^2 + Tni?p 



(27r)3 (27r)3 

g-i(£c,fci .ii -£c,fc2 ,i2 )(*^1 -*^2 ) 



fF,ii{tz2',kl)fF,i2{tz2,k2) ~ fFAi{tz2', kl) fF,i2{tz2', 1^2) 
fF,ii{tz2ikl)fF,i2{tz2,k2) - fFAi{tz2', kl) fF,i2i'tz2, 1^2) 



+ 0{tz2-tz,) 

where we set 



fjM {tzi ; ki)fF,i2 (tzi ; ^2) -fpM i^zi ; ki)fF,i2 {tzi ; ^2) 
(^zi ; ki)fF,i2 (tzi ; ^2) -fFAi {tzt ; ki)fF,i2 {tzi ; ^2) 



fF,i[t\q) = n+{t;q), fFait^l) = 1 - n-{t;q), 
fF,i{t;q) = 1 - n+{t;q), jFaii^q) = n^{t;q). 



, (227) 



(228) 
(229) 
(230) 
(231) 



It is found that the chemical potential is canceled out from the exponent in 

Eq. (mzD. 

The time evolution equation for the Bose particle is deri ved from the self- 
consistency renormalization condition for the neutral scalar fielcP2l. It is given by 



hB{tx;p) ^ -"iujpdt^ lim j dtz^dtz2Dl]j^{tj: -tz^]p)\iT}'^{tz^,tz2]p) 

+nB{tz2; p)'i^Bitzi,tz2; p) ~ nsitz^; p)iT,A{tzi,tz2] p)^ Dl'^Ritz2 ~ty]p) .(232) 

Inserting the components of the propagator (j22ip . p22p and the 1-loop quantum 
correction (|227p into Eq. (|232p . we obtain 



, EFM.ii^FM 



ii=l 12 = 1 



dt. 



d^ki d^k2 
(2^(2^ 



{27rf6^''\p^ki+k2) 



ki ■ k2 + mj 



COs{(cJb^P — EF,ki,ii + EF,k2M)i^x — ts)} 



^B,pEfM M EF,k2 ,12 
x{nB{ts;p)fF.ii{ts; ki)fFS2{is\ k2) - (1 + nB{ts;p))fF.ii{ts; ki)jF,i2{ts] k2)]. 

(233) 



This equation has the consistent statistical structure with the quantum Boltzmann 
equation. 

At 1-loop the thermal self-energy for the Dirac field, Fig. [2l^b), is given by 

d^fei d^k2 
(27r)3 (27r)3 

X r^{-^)9''{So{tz, , tz2 ; fei)r3r^(-z)/{i?o(tz2 , ^z, ; k2)T^Y\ (234) 



-^K^tz.,tz2■.p) ^ J §^§^i^-r^''Hp-k,+k2 
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Substituting the free propagators (|B.8[) and (|C.8[) and contracting the thermal index 
"7", it is found to be 

2 2 

-1} 



ii = li2 = l 

0{tzi - ^22) 



/F,ii(iz2; fcl)/s,i2(*Z2; ^2) fF,iiitz2, kl)f_B,i2itz2', ^2) 
fp Alii 7.2^1) fB,i2{iz2\k2) fF.ii{tz2;ki)fBA2(tz2;k2) 

-_/F,ii(izi; 12(^21 ; ^2) fF_,ii{tzi;ki)fBA2{izi--,k2) 

fF,iii'tziikl)fB,t2{tzi;k2) -/F,ii(izi; fcl)/B,i2(^2i; ^2), 

(235) 

where we write 

EB,q,l = ^B,q, EB,q,2 = —^B,q, (236) 

fB,i{t-,q) ^nB{t;q), /b,2(^; g) = 1 + n^lt; (?), (237) 
/B,i(t;g) = l + nB(t;g), /b,2 (*;(?) = n^lt; g). (238) 

We substitute the self-energy (|235p into Eqs. ()214p - (|215p and integrate over the 
time variables. After taking the equal time limit, we find 

lim / dtz-^dtz2S'S^p^i{tx,tzi,tz2,ty;p) 



9' 



ii=l ^2 = 1 



uJF.pEpM.ii - P • fci + ™| sin{(a;+^p - Ecmm + EB.k2.i2){tx - ts)] 

X ~~ 

2uJF,pEF,kiM'^B,k2 '^+,P ~ Ec,ki,ii + EB,k2,i2 

x{(l - n+{ts;p))fF,iAts; ki)fB,i2{ts; k2) - n+{ts]p)fF,i^{ts; fci)/s, 12(^5; ^2)}, 

(239) 

lim / dtz^dtz26Yif^^{tx,tzi,tz2-,ty]p) 

— (Pki d^k2 .3r(3), . , , ^WF,p7°+P•7-"^F 



X 



ii = l 12 = 1 

iOF,pEF,ki,ii +p ki~ nip sin{(a;„,p + Ecmm ~ -E's,fe2,J2 " 



lujF.p 



2uJF,pEF,kiAi^B,k2 ^~,p + Ec,ki,ii — EB.k2,i2 

x{n^{ts;p)fF,ii{ts;ki)fB.i2{ts; k2) - (1 - n^{ts;p))fF,n{ts; ki) f b ,i2{t s] ^2)}- 

(240) 

Due to the relationship (j229p . these equations are independent on /i. We insert 
Eqs. (|239l) - (f240l) into Eqs. (f2T4l) - (|2T5l) . and find the time evolution equations for 
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the Fermi particle and anti-particle number densities, 



(27r)3 (27r)3 
cos{(a;F,p - Efmi,ii + EBMM)iix - ts)} 



2uJF,pEF,ki,ii^B,k2 

x{(l - n+{ts;p))fF,hits; ki)fB,i2{ts; ^2) - n+{ts;p)fFMi'ts; ki) f b s] ^2)}, 

(241) 



, i^F,pEF,ki,ii + P • - 



ii=l 12 = 1 ' 

2 

COs{(wf,p + EF,ki,ii — EB,k2,i2){tx — ts)} 



2uJF,pEF,ki,iii^B,k2 

x{n-{ts;p)fFMits; ki)fBA2{ts; ^2) - (1 - n^{ts;p))jF,i^{ts] ki) f b ,i2{t s] ^2)}- 

(242) 

The last lines in Eqs. ()24ip and ()242p have the consistent forms with the collision 
term in the quantum Boltzmann equations. 



8. Conclusion 

The relativistic complex scalar and Dirac fields have been discussed at finite chem- 
ical potential in NETFD. The thermal degree of freedom is introduced through the 
Bogoliubov transformation. The fields are decomposed in terms of the Bogoliubov 
transformed operators, and Typ, with Hermitian energy eigenvalues. The time de- 
pendent Bogoliubov transformation modifies the unperturbed hat-Hamiltonian and 
induces the thermal counter term. The self-consistency renormalization condition 
has been generalized for the complex scalar and the Dirac fields. Imposing the self- 
consistency renormalization condition, we obtain the quantum Boltzmann equations 
at the leading order. A contribution beyond the quantum Boltzmann equation may 
appear from the higher order. 

In Ref. [23] the construction of a relativistic Dirac field in NETFD has not been 
accomplished by straightforward extending for the non-relativistic fields. In our 
formalism the time dependence of the particle and the anti-particle is independently 
fixed. As is mentioned in Ref. 1231 the charge conjugation invariance can be broken. 
Hence we introduce the time dependent chemical potential, /i, and successfully 
quantize. 

There are some alternative methods to derive the quantum Boltzmann equa- 
tion. In Ref. [30] the quantum Boltzmann equation for a relativistic neutral scalar 
field has been derived from the Schwinger-Dyson equation in NETFD. The thermal 
propagator in an equilibrium system has been simply extended to the one with the 
time dependent particle number density. The full propagator has been calculated 
based on the extension with the same diagram in an equilibrium system. In this 
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paper the thermal propagators are derived as an expectation value of the compos- 
ite operators constructed by two fields, (jB.l[) and (jC.l[) . This procedure based on 
the canonical quantization is more fundamental and essential in NETFD. It is now 
possible to approach high energy non-equilibrium phenomena in NETFD starting 
from the Schwinger-Dyson equation and the canonical quantization. 

There are some remaining problems. In this paper we assume that the energy 
eigenvalues for the Bogoliubov transformed free fields are independent on time vari- 
able. However, the time dependence for the eigenvalue is induced thought quantum 
corrections in an out of equilibrium system. We have to generalize the procedure 
with the time dependent energy eigenvalue at higher order. Here the thermal Bogoli- 
ubov parameters for the Dirac field are represented as spin independent quantities. 
It is necessary to extend the parameters as a spinor dependent ones for a helicity 
dependent configuration. 

A symmetry behavior is one of the most important problem in non-equilibrium 
phenomena of relativistic quark matter. It is investigated by observing an order 
parameter which is described as an expectation value of fields. Since we directly 
evaluate the expectation value of fields in NETFD, the procedure developed here 
expected to be useful for calculating the order parameter which determines the 
symmetry breaking. 

It is also interesting to apply the procedure to the early stage evolution in 
heavy ion collisions. For this purpose we should consider a spatially inhomogcneous 
state. In Refs. HOI and [H] spatially inhomogcneous states have been studied for 
the cold atom system according to the Bogoliubov-de Gennes (BdG) method in 
NETFD. The field is expanded by a complete set of wave functions under the 
spatially inhomogcneous potential. To apply the BdG method, it is necessary to 
impose the field obeys the canonical (anti-)commutation relation. It is expected that 
the procedure developed in this paper is applicable to the spatially inhomogeneous 
system. 

We hope to solve these problems and report the result in future. 
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Appendix A. Finite Chemical Potential in TFD 

We discuss a chemical potential, fj,, for complex scalar and Dirac fields in TFD. At 
fj, — a, free hat-Hamiltonian for complex scalar field is given by 




(A.l) 
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The complex scalar field is decomposed as 



C(^) = / {«p(i.)e^P-^ + {r^U^fre-^P-''} , (A.2) 

C(-) = / {S"(t.)e-P-^ + {bAt.fr.re^P-''} , (A.3) 



and 

(27r)3V 2 

' / (^/?{-«p(*^)^'^^ " + (^p(^^)''-3)"e^^-}, (A.5) 

with the following time dependence, 

a^(t,) = a^e-^'^'''-, a^(t,) = a^e"^-*^ (A.6) 

= &pe-*'^-*^ 6^(t,) = &^e*-^*^. (A.7) 

The canonical commutation relations are given by 

[C(i,a;),^f(t,y)] = mt.x)y,{t,y)]=z5^^\x-y)5^P. (A.8) 

In TFD a conserved charge, Qconsv, i-e. the zeroth component of the Noether 
current, is represented as 

Qconsv j dPx - • (A.9) 

The chemical potential, /i, can be introduced for this conserved charge. Thus the 
free hat-Hamiltonian (|A.ip is modified as 

^0 = ^^o|u=0 ~ fJ'Qconsv 



(Fx 



{ra{x)K{x) - 7f,"(x)C(a:)}] . (A.IO) 

It should be noticed that the conserved charge, Qconsv i commutes with the hat- 
Hamiltonian, iJ|p=o- It means that time evolution by the finite chemical potential is 
independent on that by H\f^=o. Then the time dependence of the ordinary operators 
for the complex scalar field (jA.2[) - (|A.5P is modified as 





= e*-f^«*-a;^(0)e 






(A.ll) 




= e*-^«*-a;^(0)e 






(A.12) 








Up c , 


(A.13) 








Up e , 


(A.14) 
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where we set 

Lo+.p = Up - ^i, (A. 15) 

LO-^p — LOp + fi. (A. 16) 

We replace the operators (|A.6|) and (|A.7p to (|A.lip - (|A.14|) . respectively. The com- 
plex scalar field (|A.2p - (jA.5P with the operators (jA.lip - (jA.14p satisfies the canonical 
commutation relations (jA.Sp . The time evolution equations of the fields (|A.2p - (jA.5P 
are obtained by the Heisenberg equations. 

In the equilibrium system the thermal Bogoliubov transformations are defined 

by 

^^{t) = B{n+{p)rPaP{t), W = aP{t)B-\n+{p)f^, (A.17) 
<W ^ S(n_(p))"^6^(i), V^it) = mt)B'\n.{p)f'^, (A.18) 

where the thermal Bogoliubov matrices, B and B~^, are given by Eqs. ([T7| and (|18p . 
In the equilibrium system the thermal Bogoliubov parameters, n±(p), for the scalar 
field represent the time independent particle and anti-particle number densities, 

where /? denotes the inverse of the temperature. 

Since the thermal Bogoliubov transformations in the equilibrium system (|A.17p 
and IjA.lSp are independent on time, the time dependence for the transformed op- 
erators, ^p(i) and rip{t), coincides with the one for the ordinary oscillators, ap{t) 
and 6°(i), 

^p{t.) = ^p"e-'"+-*^ e;(i.) = (A.21) 
'7p(ix) = ryp"e-^"--*% ^^"(t,) = jy-e^"-*^ . (A.22) 

Thus the complex scalar field is decomposed in terms of the transformed operators 
(|X2T1) and (|X22|) as 

^n^) = / {^pit^)e'P-'' + {r3fjp{t^fre-'P-^} , (A.23) 

^« ^ / (C?^7ir ^^P^^^^^'^'"' + W^f^^n'^-""} ' (A-24) 



and 



TT^ix) ^ {-^) 
n^ix) = {-^) 



I -^^^^^P^'^^'^''''' - (-3^p(^.)^)"e-^^-^} , (A.25) 
/ (^/^{-fp(*-)^"^^'' + ('?p(*-)'^-3)"e^P-^}. (A.26) 
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The complex scalar field and its canonical conjugate (jA.23|) - (IA.26|) satisfy the equal- 
time canonical commutation relations, 

= [r^{t,xUl{t,y)] ^ i5'^^\x - y)5^P . (A.27) 

Our definition of the thermal vacuum preserves the representations of the complex 
scalar field, ((X23)) - (|X26)) . in NETFD. 

In a similar manner the chemical potential can be introduced in the Dirac field. 
The conserved charge, Qconsv, is given for the Dirac field as 

d3a;V;?(x)70C(x). (A.28) 

In the equilibrium system the conserved charge introduces the chemical potential, 
fi. Thus the free hat-Hamiltonian is modified as 

Hq = Hf)\f^=o — I^Qconsv 

d^'xij^ix) • V, + m - M7°) (A.29) 
At finite chemical potential the Dirac field is decomposed as 

(A.30) 



V / "Y P s 



(A.31) 



with 



ar(^x) = a;'"e-^'^+-*% a;'"(t,) = a;■"e'"+-*^ (A.32) 

= 6;^"e'*'^"-*^ ^"{t,) = 6^'"e^'^--*^ (A.33) 

where the energy eigenvalues, uj±^p, are given by Eqs. (jA.lSp and (|A.16I) . The eigen- 
functions, and v^{p), satisfy, 

(7°cjp-7-p-mK(p) =0, (A.34) 

(-7*^^^ + 7 • p - m)i;^ (p) = 0, ( A.35) 

«^(p)(7"wp-7-p-m) =0, (A.36) 

r(p)(-7°t^j,+7-p-TO) -0, (A.37) 

and 

(P)^' (P) = -f°^P -l-p + m, ( A.38) 

S 

^ {p)v-' [p) = -f^ujp -j-p-m. ( A.39) 
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We should notice that these conditions contain no chemical potential. The Dirac 
field with the chemical potential obeys the anticommutation relation, 

{^t"(t, x)^^0^t, y)} = S^'Hx - y)S"^, (A.40) 

where we write ipl" = ■ The time evolution equations for the fields (jA.30l) and 
(|A.3ip are reproduced by the Heisenberg equations and Eqs. (|A.34p - (|A.37p . 

In the equilibrium system the thermal Bogoliubov transformation are defined 

by 

^^■"(i) = B{n+{p)rPaf{t), = -af{t)B-\n+{p)f-, (A.41) 

r;;-"(0 ^ B{n_{p)rHf{t), nrit) = bf{t)B''\n,{p))^'^, (A.42) 

where the thermal Bogoliubov matrices, B(n±) and B~^(n±), are given by 
Eqs. (|134p and (I135p . The thermal Bogoliubov parameters in the equilibrium system 
are defined with the Fermi particle and anti-particle number densities, ^ 

From the thermal Bogoliubov transformations (jA.4ip and (IA.42P we obtain the 
time dependence for the transformed operators, ^^'"(t) and rip"{t), at finite chemical 
potential, 

= C;-«e-'"+-*% ^^'"(i,) = e;^"e*'^+-*^ (A.45) 
v;'"{t.) = ?7^'"e-^"--*% %'"it,) = fjp'^e'^-^-'^. (A.46) 

With the transformed operators (jA.45p and (|A.46p we can construct the Dirac field 
at finite chemical potential as 

(A.47) 
(A.48) 

The Dirac field satisfies the canonical anticommutation relation, 

{i^l'^it, x), ^^it, y)} = S(^'\x - y)5^P. (A.49) 

In our definition of NETFD the Dirac field has the same expressions in Eqs. (jA.47P 
and (jXiS]) . 
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Appendix B. Propagator for a Free Complex Scalar Field 

In TFD the Feynman propagator for a free complex scalar field is given by the 
expectation value of the time ordered product of a composite operator constructed 
by two scalar fields. It has the 2x2 matrix form in the thermal doublet notation, 

Df{t,,ty,x-y) ^ {9\T[ra{xmy)m- (B.l) 

After the thermal Bogoliubov transformation the propagator (|B.ip reads 

DQ^{t^,ty,x - y) 

+{T3i?(n_(t,;|V,|)rr^^ 

x(^?|{r30e,_(x)P{0«,+ (y)r3p|^^){S-Hn_(t,;|V,|)rr3p'^^ 

+{T3i?(n_(i,;|V,|)fr^^ 

x(^?|{4+(y)r3rHT30«,-(a;)P|^^){S-Hn_(i,;|V,|)rr3p'^], (B.2) 

where the fields, 4''^^ and ^^j., represent the positive and the negative frequency 
parts in Eqs. ([^^ and (P5|) . respectively. 

Due to the definition of the thermal vacuum (fT9|) and \2Q\ . the thermal propa- 
gator, -Dq'^, reduces to 

Df{t^,ty, x-y) = B-\n+{t,; |V,|))"'^^Dj^^^=^(a: - \^y\)r^^ 

+ {r,Bin^it,; \S7 ^Df y''' D^,J (x - y){i3-i(n.(t,; |V,|)rr3}^^^ (B.3) 

where -Dg^j^^ (x — y) and D^^J^ {x — y) are given by 



Dl]^{x^y)= J A^^(t,-t,)-i-e— +^^(*^-*«)e^P-(=^-y) 

2u. 



DlU^ - y) 

D'oUix^y) 



other components — 0. 



{t^-ty)^ip-(x-y) 



[ ^ P o(f _f \J_pi^-.pit:o-ty) -ip-{x-y) 
J (2^)3*^^*^ ^^'2cjp^ ' 



iu-,p(t^^ty) ~ip-{x~y) 



(B.4) 
(B.5) 
(B.6) 
(B.7) 



In a homogeneous non-equilibrium system the t-representation is convenient for 
practical calculations. Performing the spatial Fourier transformation, we rewrite the 
thermal propagator in the t-representation, 

Dfit,,ty;p) - B-\n+it,;p)r-D^,^^^iU-ty;p)Bin+ity;p)r^^ 
+{T,B{n^it,;p)fr''^D^^^l^it, - ty;p){B-\n^ity;p)fT,r^^, (B.8) 
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where we define 



1 

I 

2aj„ 



-i^+,p(io:—ty) 



2(jJ„ 



Dl\{t, -ty-p)=e{ty-t 



y ^x) 



2uj„ 



Dfj,{t, - ty-p) = -e{t, - i^)-^e*"--(*^-*''), 



'0,A\''x •'u 

other components = 



2uj„ 



(B.9) 
(B.IO) 

(B.ll) 
(B.12) 



These expressions are used to evaluate Feynman diagrams for the complex scalar 
field. 



Appendix C. Propagator for a Free Dirac Field 

We evaluate the thermal propagator for a free Dirac field in NETFD. The thermal 
propagator has the 2x2 matrix form in the thermal doublet notation. It is defined 

by 



(C.l) 



Sf{t,,ty,x-y) ^ {e\T[4>^{x)i;^Jy)]\e). 
Applying the thermal Bogoliubov transformation, we obtain 
Sq^ {tx, ty, X ~ y) 

= 9{tx-ty)[B-\n+itx;\Vx\)rHo\i^j:4^)r^:-(^^^^^ 

x(^?|{r3Ve,-(x)^UV;^,+(2/)r3n^^){i3-^(n-(^^.;|Vyl))'^T3p'5" 

^6{ty-tx)[B-\n+{tx; \Vx\)r^{e\^^^%{y)i;j;4x)\9)B{n+ity; jV.DP'^ 
+{T3B(n_(t,;|V.|)m^ 

x(0|{^^,+ (y)T3rHT3V'C.-(x)rM^){S-^(n_(t,;|V,|))^r3P^], (C.2) 

where the fields, ip^j. and '>jj^±, show the positive and the negative frequency parts 
of the Dirac field, (IT5S)) and ([T571) . 

Since the transformed operators annihilate the thermal vacuum, the thermal 
propagator reads 

+{T,B{n^{U,; |V,|))^}"^^5j;^(x - y){B-\n^{ty; l^yDfr^r^^, (C.3) 
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where S'J^^^ (x — y) and S'^'^J^ {x — y) are written 

SU--y) = /(^^(^.-^.) "^^°',^^'^''"^ e-^^^(*--«)e-P-("-^), (C.4) 
Sl%{x -y) = -J ^A^y - ^.)^^^^^^^^^^^e---^(*^-»)e^P-(--?/)(C.5) 

Sl'A^-y) = J ^e{t, ^ ~ ~ "^ e--^(*--«^)e-'P-(^^^), (C.7) 
other components = 0. 
The thermal propagator (jC.3|) is expressed in the t-representation as 
^o"''(i,,t,;p) = B-\n+it,;p)r^S^^^^it,-ty;p)B{n+{ty;p)r^^ 

+{TsB{n^it,;p)fr''^S^,^2'it, - ty;p){B-'in.ity;p)fr:,r-^, (C.8) 



where 







= o{t, - 




-ty;p) 


= -Oity 


^()Ia(.^x 


-ty\p) 


- -e{ty 


^O^A i^x 


-ty;p) 





2c^ 



2ujp 



(C.9) 
(C.IO) 



_ t^fP^"+P■l--\^^-.At.-t.)^ (C.ll) 

(C.12) 



2ajp 



2u!p ' 

other components = 0. 
We use these expressions in the calculations of the 1-loop self-energy. 
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